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Density functional theory is used to study colloidal hard-rod fluids near an individual right- 
angled wedge or edge as well as near a hard wall which is periodically patterned with rectangular 
barriers. The Zwanzig model, in which the orientations of the rods are restricted to three orthogonal 
orientations but their positions can vary continuously, is analyzed by numerical minimization of the 
grand potential. Density and orientational order profiles, excess adsorptions, as well as surface 
and line tensions are determined. The calculations exhibit an enrichment [depletion] of rods lying 
parallel and close to the corner of the wedge [edge]. For the fluid near the geometrically patterned 
wall, complete wetting of the wall - isotropic liquid interface by a nematic film occurs as a two-stage 
process in which first the nematic phase fills the space between the barriers until an almost planar 
isotropic - nematic liquid interface has formed separating the higher-density nematic fluid in the 
space between the barriers from the lower-density isotropic bulk fluid. In the second stage a nematic 
film of diverging film thickness develops upon approaching bulk isotropic - nematic coexistence. 

PACS numbers: 61.30.Gd, 61.20.-p, 82.70.Dd 



I. INTRODUCTION 



There is growing interest in properties of suspensions 
of colloidal particles near structured walls because of use- 
ful applications such as selective deposition ofparticles 
and controlled growth of colloidal crystals While 
experimental HI , theoretical 0, El 0> E| j an d computer 
simulation 0>l3,Bi0| studies have been devoted to the 
understanding of the behavior of spherical colloidal par- 
ticles near geometrically structured substrates, suspen- 
sions of rodlike colloidal particles in contact with such 
substrates have not been investigated yet, despite the 
importance of rodlike colloids for both biological and ma- 
terials application. From a theoretical point of view, as 
compared with fluids consisting of spherical particles the 
study of rods is more difficult because of the additional 
orientational degrees of freedom. 

Here we study hard-rod fluids near geometrically struc- 
tured walls within the Zwanzig approximation [llj. In 
this model the allowed orientations of the rods are re- 
stricted to three mutually perpendicular orientations 
rather than a continuous range of orientations in space 
(see Fig. ^| ; the positions of the rod centers are continu- 
ous variables. The advantage of this model is that the dif- 
ficult determination of spatially inhomogeneous density 
and orientational order profiles becomes feasible allow- 
ing one to study various aspects of hard-rod fluids near 
structured walls in detail. On the basis of recent theoret- 
ical studies on fluids of hard rods near planar hard walls 
UH IT3l Il4| , we expect to find results which remain qual- 
itatively correct even in the absence of the restriction to 
discrete orientational directions. In Sec. II we describe 
the density functional theory which is used to analyze a 
hard-rod fluid in contact with an individual right-angled 



wedge or edge (Sec. Ill) or with a periodically patterned 
wall (Sec. IV). 



II. MODEL AND DENSITY FUNCTIONAL 
THEORY 



We consider an inhomogeneous fluid consisting of hard 
rods of length L and diameter D. The number density 
of the centers of mass of the rods at point r with orien- 
tation uj = (0, (/)) of the normal along their main axis of 
symmetry is denoted as p(r,o>). The equilibrium density 
profile of the inhomogeneous liquid under the influence 
of an external field V (r, oj) minimizes the grand potential 
functional 

Q\p(r,w)] = J d 3 rdup(r,u) [k B T (ln[47rA 3 p(r, lu)} - l) 
- -Li + V{r,u)]+F ex [p{r,uj)], 

(1) 

where A is the thermal de Broglie wavelength and /1 is 
the chemical potential. Within the Onsager second virial 
approximation the free energy functional F ex [p{v,uj)\ in 
excess of the ideal gas contribution reads [I3 



F EX [p(r,uj)} = -~y~ J d 3 r 1 dL0id 3 r 2 duj2 piri,^) 



xp(r 2 ,o;2)/(ri,r2,wi,L i ;2) 



(2) 



where /(iq, r 2 , UJ\, UJ2) is the Mayer function of the inter- 
action potential between two rods. The Mayer function 
equals —1 if the rods overlap and is zero otherwise. 

In the present application of density functional the- 
ory we concentrate on the ordering effects induced by 
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FIG. 1: The system under consideration consists of a fluid of 
thin hard rods of length L and thickness D <C L in contact 
with a hard substrate (grey). The surface z s (x) of the sub- 
strate wall exhibits a periodic pattern with period p consisting 
of rectangular blocks of width w and height h. The density 
profiles for the centers of the rods with an orientation of the 
normal along their main axis of symmetry parallel to the or- 
thogonal unit vectors of the reference frame are denoted as 
p x (x, z), p y (x, z), and p z (x, z). The system is spatially invari- 
ant in {/-direction. In accordance with the discreteness of the 
orientational degrees of freedom p z (x,z < z s (x) + L/2) = 0, 
p x (x + np - (w + L)/2 < x < xo + np + (w + L)/2, z < 
h) = with n € Z, and p X]V (x,z < z s (x)) = where 
z s (xo + np — w/2 < x < xo + np + w/2) = h and z a (x) = 
otherwise. The value of xo is arbitrary. 

surfaces geometrically structured such that the result- 
ing p(r, oj) depends on two spatial coordinates. For the 
model system displayed in Fig.^ apart from the possibil- 
ity of surface freezing at high densities, non-uniformities 
of the density occur only in the x-z plane, so that 
p(r, oj) = p(x, z, 9, </>). Minimization of f2 with respect to 
p{x, z, 9, (f) leads to the following Euler-Langrange equa- 
tion: 

k B T\n[4TrA 3 p(x, z, 9, 4>)] = p - V{x, z, 9, <f>) 

+ J dx\dz\ J dcf>i J d9\ sin 6*1 p(x\, z%, 9\, 4>\) 
o o 

xf(x,z,6,<f),x 1 ,zi,6x,(pi)kBT. (3) 

This equation can be solved numerically for a given chem- 
ical potential p and a given external field V(x, z, 9, <fi). 
For computational purposes, the density profile has to be 
specified on a sufficiently fine four-dimensional (x, z,9, <p) 
grid. In order to reduce this c omp utational effort we use 
the Zwanzig model for rods pT|. Within the Zwanzig 
model the rods are represented by rectangular blocks of 



size L x D x D. The positions of the center of mass 
vary continuously, while the allowed orientations of the 
normal of each rod are restricted to directions parallel 
to the x, y, and z axis (see Fig. QJ. Using the notation 
a x (x, z) — a(x, z,9 = 7r/2, (f> = 0), a y (x, z) = a(x, z,6 = 
7r/2, 4> = tt/2), and a z (x, z) = a(x, z,9 = 0,<fi = 0) with 
a = p, V, the Euler-Langrange equations for a fluid con- 
sisting of thin Zwanzig rods {D/L — > 0) can be written 
as 

lnlA 3 ^^,^)] = {k B T)- l {p~V x {x,z)) 

px+L/2 pz+L/2 

-2D / dxi / dz\ p z (xi, zi) 

Jx-L/2 Jz~L/2 
rx+L/2 

-2DL / dx 1 p y (x 1 ,z) , (4) 

Jx-L/2 

In[A 3 p„(a;,«)] = {k B T)-\n - V y {x, z)) 

rz+L/2 

—2DL j dz\p z {x,Z\) 

Jz-L/2 
rx+L/2 

-2DL dx 1 p x (x 1 ,z), (5) 

Jx-L/2 

and 

]n[A 3 p z (x,z)} = {k B T)- 1 {ii-V z {x,z)) 

px+L/2 pz+L/2 

-2D j dxi j dzi p x (xi, zi) 

Jx-L/2 Jz-L/2 
i-z+L/2 

-2DL / dz 1 p y {x, z\) , (6) 

Jz-L/2 

which allow for a straightforward iterative numerical 
computation on a two-dimensional (x, z) grid. We note 
that the meaning of the terms on the right side of 
Eqs. Q)-© can easily be inferred from considering the 
various orientations of the rods (see Fig. ^) . It is conve- 
nient to introduce the variable p* = p— k B T hi(A 3 / L 2 D) . 
In the following sections numerical data are given in 
terms of p* and we drop the star in order to avoid a 
clumsy notation. 

The bulk phase behavior of this model was studied a 
long time ago by Zwanzig pT|. who found a first-order 
isotropic - nematic phase transition similar to Onsager's 
result for freely rotating rods 0] . Recently, the Zwanzig 
model has been used to investigate the phase behavior of 
monodisperse |l2l ll3T ] and binary [l4T | rod fluids near a 
single planar hard wall and confined in slit pore. These 
calculations yield a wall-induced continuous surface tran- 
sition from uniaxial to biaxial symmetry. Complete wet- 
ting of the wall - isotropic liquid interface by a biaxial ne- 
matic film has been found. For the fluids confined by two 
parallel hard walls, at large slit widths a first-order capil- 
lary nematization transition occurs, which terminates in 
a capillary critical point upon decreasing the slit width. 
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FIG. 2: Cross-section of an individual right-angled wedge and 
edge, respectively. The outer part of the shaded region forms 
an edge whereas the inner part forms a wedge. The geometri- 
cal structures are translationally invariant in the y-direction 
perpendicular to the plane of the figure. In Figs. 0|and|^|den- 
sity profiles of thin hard rods of length L are shown along the 
paths indicated by the solid, dashed, and dotted lines which 
are parallel to the outer and inner surface, respectively. The 
paths are parameterized by s such that for each line s = 
(•) indicates the position of its corner at x = z. p y {x,z) is 
symmetric around s — and p x for s ^ equals p z for s «s 0. 
Both for the wedge and edge the position of the corner is 
associated with (x,z) = (0,0); in this figure both cases are 
superimposed. 



III. HARD-ROD FLUID NEAR A 
RIGHT-ANGLED WEDGE AND EDGE 

Before studying the hard-rod fluid near the geometri- 
cally structured surface shown in Fig.^it is instructive to 
analyze first the fluid around an individual right-angled 
wedge or edge (see Fig. These simple geometrical 
structures constitute the building blocks of the structured 
surface displayed in Fig. ^ An analysis of the size depen- 
dence leads to the following decomposition of the grand 
canonical potential functional of the fluid which in its 
bulk is taken to be in the isotropic phase: 

U[p{x, z)] = H y [H x H z io b + (H x + H z )j wI + 7-/(0!)] , (7) 

where u>b is the bulk grand canonical potential density, 
j w l is the wall - isotropic liquid surface tension at a pla- 
nar wall and 77(a) is the line tension of the isotropic 
liquid with a = ir/2 for the wedge and a = 37r/2 for 
the edge. The extension Hp of the system in direc- 
tion (3 — x,y,z is defined as the length available to the 
rim of the particles at closest approach to the boundary. 
We restrict our analysis to chemical potentials ji smaller 
than the chemical potential ^ IN) /{k B T) = 0.8087227 
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FIG. 3: (a) The wall - isotropic liquid surface tension at a 
planar wall j w i (solid line) together with the line tension at a 
right-angled wedge 77(71-/2) (dashed line) and edge 77(371-/2) 
(dotted line) of a fluid consisting of thin rods of length L and 
diameter D (D/L — > 0) as a function of the chemical poten- 
tial. The vertical line marks the location of the chemical po- 
tential fi (IN) /(k B T) = 0.8087227 at bulk isotropic - nematic 
coexistence. The wall - isotropic liquid interface is completely 
wetted by a nematic film, i.e., j w i — j w n + Jin — 0.0498 
k B T/(LD) at /i (/Ar) , where j w n is the wall - nematic liq- 
uid surface tension and jin is the isotropic - nematic in- 
terfacial tension E3, HE E3 At p {IN) , Ti(3n/2) is slightly 
larger than j w i L and 77 (tt/2) attains a finite value which 
becomes visible only at higher resolutions. 77 (tt/2) changes 
sign at /i/(fcsT) = 0.27. (b) The total particle number den- 
sity pb = p(x, z —* 00) [see Eq. JHJ] of the homogeneous and 
isotropic bulk fluid as a function of the chemical potential: 
p/(k B T) = ln(p b DL 2 ) -ln3 + 4p b L>L 2 /3- At bulk isotropic - 
nematic coexistence (marked by the vertical line) the density 
of the isotropic phase is given by p b DL 2 = 1.25822486. This 
figure allows one to translate values for p, into pb, which is 
experimentally accessible, and vice versa. 



at bulk isotropic - nematic coexistence. Figure dis- 
plays the surface tension j w i as well as the line tensions 
77(77/2) and 77(37r/2) as a function of the chemical po- 
tential. The steric interaction between the particles in- 
creases the surface tension with increasing chemical po- 
tential. On the other hand, the onset of the surface- 
induced nematic ordering of the particles leads to a de- 




FIG. 4: The density profiles p x (x,z), p y (x,z), and p z (x,z) for 
thin hard rods of length L and diameter D (D/L — > 0) near 
a right-angled wedge evaluated along the paths specified in 
Fig.|5][using the same line code, i.e., solid, dashed, and dotted 
lines]. The chemical potential is fi = 0. The solid curve in 
(a) [(c)] represents for s > [s < 0] the density profile along 
the line of closest contact, i.e., rods oriented parallel to the x 
[z] axis touch the wall with the rim. p y is symmetric around 
s — and p x for s ^ equals p z for s «j 0. 



FIG. 5: The density profiles p x (x,z), p y (x,z), and p z (x,z) for 
thin hard rods of length L and diameter D (D/L — * 0) near a 
right-angled edge evaluated along the paths specified in Fig. [21 
[using the same line code, i.e., solid, dashed, and dotted lines]. 
The chemical potential is p = 0. The solid curve in (a) [(c)] 
represents for s > L/2 [s < —L/2] the density profile along 
the line of closest contact, i.e., rods oriented parallel to the x 
[z] axis touch the wall with the rim. p y is symmetric around 
s — and p x for s ^ equals p z for s «j 0. 



crease of the surface tension for larger chemical poten- 
tials. In the limit of large negative chemical potentials, 
i.e., for non-interacting particles, the wall - isotropic liq- 
uid surface tension as well as the line tensions vanish. 
The line tension tj(tt/2) for the fluid near a right-angled 
wedge exhibits a change of sign with increasing chemi- 
cal potential while tj(3w/2) > for all values of pi. We 
note that by construction Ti(jt) = 0. As a function of a 
the line tension 17(a) corresponds to the work done per 
unit length, against the fluid, to change the dihedral an- 
gle from 7r to some value a |16l Il7| . The corresponding 
solvation torque is i/(cv) = —drj(a)/da. 

Figures 21 and show the density profiles at the right- 
angled wedge and edge for a chemical potential /1 = 0. 
The profiles are evaluated along lines parallel to the con- 
fining walls (see Fig.|2J), using a linear parametrization s 
such that s = corresponds to the position of the corner. 
For s — * ±00 the density profiles at a planar hard wall are 
recovered. As is apparent from Fig. 0] the density profile 



p y {x, z) exhibits a maximum at s — for the right-angled 
wedge. This maximum increases upon approaching the 
corner of the wedge. For a right-angled edge a significant 
depletion of rods lying parallel to the y axis is found near 
the corner (see Fig. |SJ) ■ The sharp features of the den- 
sity profiles in the presence of the right-angled edge are 
caused by the vanishing of p x (x > —L/2,z > 0) and 
Pz{x > 0, z > —L/2) which is due to the presence of the 
impenetrably hard walls. Such cusps and discontinuities 
- although less pronounced - have already been found 
for the density profiles near a planar hard wall using the 
Zwanzig model [l2L IT3L Il4| , while the interfacial profiles 
between demixcd fluid phases, as calculated for the same 
model, exhibit neither cusps nor discontinuities |l8j| . We 
do not expect to observe discontinuities in the density 
profiles for freely rotating rods near an edge, similar to 
our findings for hard rods near a planar hard wall [l9| 
(see discussion below). Although the calculated density 
profiles presented in Figs. 01 and can only be consid- 
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ered to be of qualitative significance we expect that the 
main features, namely the enrichment [depletion] of rods 
lying parallel and close to the corner of a wedge [edge] 
as well as the asymmetry of the density profiles p x (x, z) 
and p z (x, z) with respect to the line s = 0, remain valid 
for freely rotating rods. 

In order to understand the structure of the density pro- 
files displayed in Fig. 0] and [SI it is instructive to apply 
the idea of entropically driven forces (see, e.g., Ref. |2jj) 
to the Zwanzig model. We consider rods of a given ori- 
entation along the x, y, and z axis as belonging to one of 
three " species" . Such a three-component fluid maximizes 
its entropy by maximizing the volume accessible per rod. 
Although there exist only steric repulsions between pairs 
of particles, maximizing the entropy in the fluid mixture 
can lead to an effective entropic attraction between the 
rods and the walls. FigurcEl(b) demonstrates that, when 
a rod of a given species approaches a planar wall (rep- 
resented in grey), the total volume available to rods of 
the other species increases. This increases the total en- 
tropy of the mixture by an amount proportional to the 
size of the excluded-volume overlap region (represented 
in black) multiplied by the pressure. For a rod lying 
close and parallel to a right-angled wedge [edge] the cor- 
responding excluded-volume overlap region is increased 
[decreased] (see Figs. El (c) and (d), respectively) leading 
to an enrichment [depletion] of such rods close to the cor- 
ner of a wedge [edge] . The results may be interpreted in 
terms of a repulsive barrier of an effective potential re- 
pelling a rod, which is oriented parallel to the corner of an 
edge, approaching the edge from the side and practically 
preventing it from passing around the corner. On the 
other hand the effective potential acting on a rod which 
is oriented parallel to the corner of a wedge is " pushing" 
it into the corner. For a detailed analysis of these mecha- 
nisms acting on mixtures of hard spheres near edges and 
wedges see Ref. [j| and in particular Figs. 5-7 therein. 

The simple illustration in Fig. (d) is also helpful for 
understanding the aforementioned discontinuities in the 
density profiles near the right-angled edge (see Fig. 0. 
When a thin rod (D/L — > 0), which is oriented paral- 
lel to the edge, approaches the edge from the side, the 
excluded-volume overlap drops abruptly to zero before 
the rod is passing around the corner. This causes the dis- 
continuities in the density profiles along the paths spec- 
ified in Fig. |21 For freely rotating rods the correspond- 
ing excluded-volume overlap decreases smoothly to zero 
because of the huge number of differently oriented rods 
acting on the rod which is oriented parallel to the edge. 

Finally, we briefly discuss the phase behavior of hard 
rods confined in a hard pore of square cross-section. The 
immediate consequence of the pore is that rods oriented 
perpendicular to the confining walls cannot approach 
closer than a center-of-mass distance L/2. There is a pro- 
nounced increase of the density of rods orientated parallel 
to the main axis of the pore in the corners of the pore be- 
cause of the aforementioned effective entropic attraction. 
For sufficiently large cross-sections of the pore, we ob- 




FIG. 6: Illustrations of the effective entropic interactions of 
hard rods with surfaces, (a) Schematic side view of rectangu- 
lar rods of length L and thickness D with orientations along 
the x, y, and z axis, where the t/-axis is perpendicular to 
the plane of the figure, (b) Due to the steric interaction, the 
centers of mass of rods with ^-direction are excluded from 
the hatched regions surrounding the rods with y-direction 
(quadratic cross-section) and the walls (thick grey lines) . Here 
one rod with x-direction touches a rod with y-direction (left) 
while another rod with x-direction fully touches the wall 
which is oriented parallel to the z-axis (right). The rods with 
x-direction are viewed as depletion agents leading to an ef- 
fective interaction between the rod with y-direction and the 
walls. When rods with orientations along the y axis are suf- 
ficiently separated from each other and from the walls, the 
volume accessible to the rods with ^-direction is the total vol- 
ume of the container minus the volume of the hatched regions. 
However, when a rod oriented parallel to the i/-axis is close 
to a wall the volume accessible to the rods with z-direction 
increases by the excluded-volume overlap region marked in 
black. The corresponding increase in entropy induces an effec- 
tive attractive force between rods with an orientation parallel 
to the t/-axis and the wall, (c) In a corner of a right-angled 
wedge (thick grey lines), the overlap volume (black region) is 
larger than the one on a flat wall leading to an enrichment of 
rods lying parallel and close to the corner of the wedge [see 
Fig. 0](b)]. (d) Close to a right-angled edge (thick grey lines), 
the overlap volume is smaller than the one on a planar wall. 
Therefore the density of rods lying parallel and close to the 
edge is smaller than the density near a planar hard wall at 
the same chemical potential [see Fig. El(b)]. Similar consid- 
erations hold for the rods with z-direction acting as depletion 
agents on the rod with t/-direction. The rods with y-direction 
are exposed to the superposition of both effective interactions. 



serve coexistence between an isotropic phase and a cap- 
illary condensed nematic phase. The density profile of 
the capillary condensed nematic phase is characterized 
by a nematic phase throughout the pore, whereas the 
density profile of the coexisting phase decays toward an 
isotropic phase in the middle of the pore. For small pore 
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cross-sections a sharp capillary nematization transition 
no longer occurs and is replaced by a steep but continu- 
ous filling upon increasing the chemical potential. For the 
same fluid confined in a slit pore the confinement effects 
are weaker. Thus, in the slit pore we observe capillary 
nematization at a higher chemical potential correspond- 
ing to a higher particle number density of the bulk fluid. 
However, the spatially averaged particle number density 
of the coexisting inhomogeneous isotropic phase in the 
slit pore is smaller than the corresponding one in the 
pore of square cross-section. 

IV. HARD-ROD FLUID IN CONTACT WITH A 
PERIODICALLY STRUCTURED HARD WALL 

We now turn our attention to the properties of the 
hard-rod fluid in contact with the hard wall shown in 
Fig. n The surface of the wall is periodically patterned 
with rectangular hard barriers of width w and height h, 
where the periodicity is denoted by p. We focus on the 
numerically determined orientationally averaged number 
density profile 

p(x, z) = p x (x, z) + p y (x, z) + p z (x, z) (8) 

and the excess adsorption Y defined as 

T = J dx dz [p(x, z) — pb] , (9) 

where pb — p(x, z — > oc) is the total particle num- 
ber density of the homogeneous bulk fluid. The volume 
V = J dxdydz of the system is defined as the total vol- 
ume of the container, i.e., the left boundary of the system 
displayed in Figgis taken to be the surface z s (x) of the 
substrate wall which implies that the trenches between 
the barriers contribute to V. Figure displays T for var- 
ious values of the barrier height h and two values of the 
barrier width w at a fixed periodicity p of the surface pat- 
tern. For non-interacting rods (p — > — oo), the calculated 
excess adsorption reveals a slight depletion close to the 
surface (r < 0) because there is less space available to the 
rods in the presence of the impenetrably hard walls. For 
the same reason this depletion becomes more pronounced 
with increasing height of the barriers (i.e., increasing the 
actual exposed solid area) . Upon increasing the chemical 
potential, the excess adsorption increases and exhibits 
a change of sign because of the aforementioned entropic 
attraction between the rods and the surface [see Fig. H5| . 
For small barrier heights, T increases smoothly upon in- 
creasing the chemical potential, while a pronounced vari- 
ation of the excess adsorption is found for large barrier 
heights at a chemical potential smaller than the chemical 
potential p 1 -- 1 ^ at bulk isotropic - nematic coexistence. 
Moreover, the calculation renders T to diverge logarith- 
mically as p — » pS IN h Near p( IN ) the excess adsorption 
can be fitted byT — Ai-A 2 ln[(/i( /JV ) - p) /{k B T)\ , with 
fit parameters Ai and A 2 , where A 2 > turns out to be 
independent of the surface pattern. 
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FIG. 7: The excess adsorption V [see Eq. (0] of a fluid consist- 
ing of thin hard rods of length L and diameter D (D/L — > 0) 
near a geometrically structured wall as shown in Fig. Q for 
various barrier heights: h — (solid curves); h = 2L (dashed 
curves); h — 8L (dash-dotted curves); h — 12 L (dotted 
curves) . The periodicity is p = 5 L and the barrier width 
is w = 3 L in (a) and w = L in (b) . The vertical lines mark 
the value of the chemical potential p (IN) /{k B T) = 0.8087227 
at bulk isotropic - nematic coexistence. In all cases T diverges 
logarithmically as p — > p^ IN K Figure |H| exhibits density pro- 
files for the system with h = 8 L in (b) for the three chemi- 
cal potentials p — 0.75, 0.77, 0.808 ksT marked by arrows. 
The chemical potential p ~ 0.757 feT corresponding to the 
pronounced variation (no jump but steep increase) of V for 
h — 12 L in (b) turns out to agree with the chemical poten- 
tial at the occurrence of the first-order capillary nematization 
transition of the same fluid confined in a corresponding slit 
pore of width p — w — 4L. Neither the curves in (a) nor in 
(b) intersect at a single point. 

The logarithmic divergence of T is consistent with com- 
plete wetting of the wall - isotropic fluid interface by a 
nematic film in the absence of algebraically interaction 
potentials |21j . A similar behavior of the excess adsorp- 
tion close to /i^^- 1 has been found for the same fluid near 
a planar hard wall dill El. 

To understand the origin of the calculated excess ad- 
sorptions, it is instructive to study the variation of the 
density profiles with increasing chemical potential. The 
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orientationally averaged density profiles p(x, z) shown in 
Fig. [S] demonstrate that the wetting of the non-planar 
wall - isotropic liquid interface by a higher-density ne- 
matic film occurs as a two-stage process where first the 
nematic phase fills the space between the barriers un- 
til an almost planar isotropic nematic liquid interface 
has formed separating the higher-density nematic fluid 
in the space between the barriers from the lower-density 
isotropic bulk fluid. In the second stage a nematic film 
of diverging film thickness develops upon approaching 
the chemical potential at bulk isotropic - nematic coexis- 
tence. In the presence of the patterned wall, the director 
(the average orientation of the rods) of the nematic phase 
is parallel to the y-axis because of the aforementioned 
effective entropic attraction between rods oriented paral- 
lel to the right-angled wedges of the barriers [see Fig. 
(c)]. Upon approaching the chemical potential at bulk 
isotropic - nematic coexistence, i.e., p > 0.757 fcgT in 
Fig. [7| (with this value of the kink position being largely 
independent of h), the calculated density profiles at the 
isotropic - nematic interface become virtually indistin- 
guishable from the free isotropic - nematic interface be- 
tween coexisting bulk phases, as expected for the case of 
complete wetting. In this limit this holds irrespective of 
the actual values of the width w, height h, and period- 
icity p. However, the thickness of the emerging nematic 
film depends on the aforementioned model parameters 
and the chemical potential. 

A very thin nematic liquid layer - corresponding to 
a large undersaturation - follows the substrate pattern, 
whereas a sufficiently thick layer is essentially flat (see 
Fig. |SJ). In the following, we will be exclusively con- 
cerned with thick nematic films, allowing us to define 
an a;- independent film thickness t = Ix/(j,|,|, rfr) 
defined as the distance z m between the midpoint of the 
density profile for the planar isotropic - nematic interface 
and the substrate surface z = 0; / is a scaling function 
appearing on the basis of dimensional analysis. The min- 
imal film thickness for which a planar isotropic - nematic 
liquid interface is still possible is t > h. We consider the 
situation as illustrated in Fig. [§| (a), where for a given 
chemical potential p, a nematic (N) film of thickness to 
intrudes between a planar hard wall and an isotropic (I) 
bulk fluid. At the same chemical potential the film thick- 
ness t of the fluid in contact with a geometrically pat- 
terned wall is larger than the height of the barriers h and 
smaller than to + h (see Fig. 03(b)). For p = w it follows 
that t = to + h as expected on physical grounds. It is 
worthwhile to mention that /(■£, j;, jr^f ) as function 
of w/L for given h/L, p/L, and fi/(kgT) exhibits a dis- 
continuity upon approaching w — ► because the ratio of 
the actual substrate area (including the side planes of the 
rectangular blocks) per period p over the one projected 
onto the x — y plane drops abruptly from (p + 2/i)/p to 1 
for w = 0. Hence, due to geometric constraints the film 
thickness in the presence of infinitely thin barriers differs 
from the one in the absence of the barriers. Moreover, 
we find the following properties, which are schematically 




z/L z/L 



FIG. 8: Orientationally averaged total density profile p(x, z) 
[see Eq. 0] of thin rods of length L and diameter D (D / L — > 
0) in contact with a hard wall for three chemical potentials 
p = 0.75, 0.77, 0.808 k B T, where p {IN) /{k B T) = 0.8087227 
is the chemical potential at bulk isotropic - nematic coex- 
istence. The corresponding contour plots are shown on the 
right. The sharp structures in the contour lines are caused by 
the discreteness of the orientational degrees of freedom. The 
wall at 2 = is patterned with rectangular (parallel) barriers 
of width w = L with xo = 2.5 L, height h — 8 L, and periodic- 
ity p = 5 L (see Fig. . The corresponding excess adsorption 
is represented by the dash-dotted line in Fig. [7| (b) , and the 
three chemical potentials are marked by arrows in Fig. Q (b) . 
For p = 0.808 k B T the interface between the higher-density 
nematic film on the wall and the lower-density isotropic bulk 
fluid [z oo) resembles closely the interface between the 
free isotropic - nematic interface between coexisting bulk 

( I ) 9 

phases, with coexisting densities p b DL = 1.25822486 and 
p[ N) DL 2 = 1.91544377. At p ~ 0.757 k B T the trenches un- 
dergo a nematization filling which is a smooth but steep vari- 
ation of the density distribution. 

visualized in Figs. 03 (b) and (c): 

.,h p w a . ,h w a . .. 

t=ix %k'^ =£x ^?^ (10) 

and for fixed w/p and /i/(/csT) 

t(h + hi) = hi + t(h) , h>L. (11) 
Equation (|l(Jfl states that upon varying the barrier width 
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FIG. 9: Schematic side view of a planar wall and two geomet- 
rically patterned walls (dark grey regions). The y-direction 
is perpendicular to the plane of the figure (see Fig. 0. (a) 
The planar wall - isotropic (I) liquid interface is wetted by a 
nematic (N) film of thickness to- (b) At the same chemical 
potential the film thickness t on the geometrically patterned 
wall is larger than h and smaller than to + h, where i = to + h 
for w = p. (c) Illustration of Eqs. 1101 and 1111 : w and p are 
half as large as in (b) and h is increased by hi. 



w and the periodicity p such that the ratio (p — w)/w of 
the substrate area at the bottom z = over the substrate 
area at the top z — his kept fixed, the film thickness does 
not change for a given barrier height h, chemical poten- 
tial and L. This result is reminiscent of the Cassie 
equation |2^ . Cassie considered a simple liquid in con- 
tact with a smooth but chemically striped surface with 
periodicity p such that on stripes of width w one has a 
contact angle d\ and in between $2- The apparent aver- 
age contact angle is given by 



cos $ 



app — — COS$i 
p 



1 COS $2 

P , 



(12) 



Since i9 a pp, $i, and $2 are determined uniquely by t app , 
ti, and £2 via the corresponding effective interface po- 
tentials (see Eq. (4.56) in Ref. |2l|), Eq. {TJJ) states that 
the apparent film thickness t app depends only on the ra- 
tio w/p, which is the analogue of Eq. (|10fl . Whereas 
Eq. ifTUI) is valid for all values of w/L, h/L, and p/L, 
for molecular-scale surface patterns with barrier heights 
h < L we find deviations from Eq. Qllfl. 

The wetting of the non-planar wall - isotropic liquid 
interface by a nematic film is driven by the steric in- 
teraction of the rods with the solid substrate which is 
mediated by the fluid occupying the space between the 
barriers. As the barrier height h increases the interac- 
tion between rods located at z > h and the part of wall 
which is located at z = weakens, such that for large 
height h and small width w, wetting is dominated by the 
interaction with the fluid in the space between the bar- 
riers and not by that with the solid substrate at z = 0. 
However, one has to take into account that the rods in 
the trenches between the barriers interact not only with 
the part of the substrate wall which is located at z = 
but also with the side planes of the rectangular blocks 
which are parallel to the z — y plane. 



V. SUMMARY 

We have studied hard-rod fluids near geometrically 
structured walls using Zwanzig's model of square paral- 
lelepipeds with only three allowed orientations (Fig. 
Within the framework of a density functional theory, the 
grand potential functional is minimized numerically and 
density profiles, excess adsorptions as well as surface and 
line tensions are determined leading to the following main 
results: 

(1) The line tension for the isotropic fluid in contact 
with a right-angled wedge (see Fig. |2J) exhibits a change 
of sign with increasing chemical potential while the line 
tension for the fluid in contact with a right-angled edge 
as well as the wall - isotropic fluid surface tension at a 
planar hard wall are positive (Fig.[3J). 

(2) Figures 0] and [5] demonstrate an enrichment [de- 
pletion] of rods lying parallel and close to the corner 
of a right-angled wedge [edge]. On the basis of effec- 
tive entropic forces between the rods and the walls (see 
Fig. [fjj) , the results may be interpreted in terms of a re- 
pulsive barrier of an effective potential repelling a rod, 
which is oriented parallel to the corner of an edge, and 
approaches the edge sidewise, and practically preventing 
it from passing around the corner. The effective potential 
acting on a rod which is oriented parallel to the corner 
of a wedge is larger than the one close to a planar wall. 
Building on the effects demonstrated in Figs.|3J|SJ and El 
it seems possible to devise structures that create local- 
ized and directional entropic force fields for both natural 
and synthetic rodlike colloids. 

(3) Coexistence between an isotropic and a capillary 
condensed nematic phase is observed for the fluid con- 
fined in a hard pore of square cross-section, provided the 
cross-section is sufficiently large. The density profile of 
the capillary condensed nematic phase is characterized 
by a nematic phase throughout the pore, whereas the 
density profile of the coexisting phase decays towards an 
isotropic phase in the middle of the pore. For the same 
fluid confined in a slit pore the confinement effects are 
weaker, i.e., in the slit pore one observes capillary nema- 
tization only at a higher chemical potential. 

(4) From the calculated excess adsorptions (Fig.|JJ and 
density profiles (Fig. ISJ of a fluid consisting of hard rods 
near the geometrically structured wall shown in Fig. ^ 
we conclude that complete wetting of the non-planar wall 
- isotropic liquid interface by a nematic film occurs as a 
two-stage process. In the first stage the nematic phase 
fills the space between the barriers until an almost planar 
isotropic - nematic liquid interface has formed separating 
the higher-density nematic fluid in the trenches between 
the barriers from the lower-density isotropic bulk fluid. 
In the second stage a nematic film of diverging film thick- 
ness develops upon approaching the chemical potential 
at bulk isotropic - nematic coexistence. The film thick- 
ness, defined as the distance between the midpoint of the 
density profile for the almost planar isotropic - nematic 
interface and the substrate bottom at z — 0, is larger for 
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the fluid near the geometrically structured wall than the 
one for the fluid near a planar wall at the same chemical 
potential (Fig. EJl- 

Finally, we note that phenomena which emerge from 
the contact of a rod fluid which is in its bulk in the nc- 



matic phase are also interesting because of the possibility 
to deliver external lateral structures deep into the bulk 
of the adjacent fluid which offers a convenient means to 
image patterned surfaces. Density functional theory will 
allow one to study also such a system. 
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